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VT -symmetric Deformations of the Korteweg-de Vries 
Equation 



Abstract: We propose a new family of complex PT-symmetric extensions of the 
Korteweg-de Vries equation. The deformed equations can be associated to a sequence 
of non-Hermitian Hamiltonians. The first charges related to the conservation of mass, 
momentum and energy are constructed. We investigate solitary wave solutions of the 
equation of motion for various boundary conditions. 



1. Introduction 

PT-symmetry has served as a very fruitful guiding principle to identify potentially interest- 
ing non-Hermitian Hamiltonians, which may constitute physically relevant non-dissipative 
systems. The interest in these type of configurations has started with a numerical obser- 
vation made in Q, where it was found that the Hamiltonian 



possess a real, positive and discrete eigenvalue spectrum for integers N > 1 with coupling 
constant g E M + , despite it being non-Hermitian H ^ and unbounded from below, for 
N = An — 1 with n £ N. The virtue of PT-symmetry results from the fact that whenever 
the Hamiltonian and the wavefunctions are left invariant under a PT-transformation the 
eigenvalues are guaranteed to be real. However, the anti-linear nature of the PT-operator 
is responsible for the fact that such a guarantee can not be provided by the PT-symmetry 
of the Hamiltonian alone [J. Unlike as for linear operators, for the PT-operator its 
two dimensional representation can be realized, in which case one speaks of broken PT- 
symmetry. One is then in a situation in which the corresponding wavefunctions are not 
PT-symmetric and the eigenvalues occur in complex conjugate pairs. Nonetheless, even 
though PT-symmetry of the Hamiltonian can not guarantee the reality of the spectrum, it 
pre-selects a subclass of promising non-dissipative systems. For recent results and a review 
see for instance || |6| . 
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A couple of month ago Bender, Brody, Chen and Furlan |?| have applied the above 
principle to identify interesting extensions of the Korteweg-de Vries (KdV) equation || 

u t + uu x + u xxx = 0. (1.2) 

The scaling properties of this equation for x — ► ax, t — > /3t,u —>■ are well known, e.g. Q 

2 

—Ut H uu x + -^u xxx = 0, (1.3) 

pa a 6 

and it has been remarked already at least thirty years ago that the KdV equation remains 
invariant under a PT-transformation t — > —t,x — ► — x, see for instance p. 414 in Q. This 
is of course just the particular case a = (3 = —7 = —1. However, this property has only 
been exploited in the above mentioned spirit in 0, where the KdV equation has been 
extended to the complex domain in a PT-symmetric manner 

u t - iu{iu x f + = eeR. (1.4) 



One may think of equation ( |1.4| ) as being obtained from ( |1.2| ) by a scale invariant defor- 
mation 

u x -> -q{qu x f £Gt, (1.5) 

of the second term. When the deformation parameters scale as q — > a/~fq, q — ► ~f/aq, 



equation (1.4) has the same behaviour under scaling as ( |1.3| ) for all values of e. The special 
case q = q = i yields a "PT-symmetric expression for a = (3 = —7 = — 1. Intriguingly, the 
equation (|l.4[) were found to possess interesting solitary wave solutions and two conserved 
charges were also constructed. 

2. A new PT-symmetric deformation of the KdV equation 

It should be mentioned that complex extensions of the KdV equation have been studied 



before, see e.g. ]10|, |ll|, |12[] and in passing even some special cases of equation (|L4j) have 



been dealt with for instance in [13]. However, only few properties have been studies for 
the latter and PT-symmetry has not been adopted as a guiding principle. Motivated by 
the interesting findings in fm and the usefulness of PT-symmetric complex deformations in 
other contexts, see e.g. [||, ||, ||, we extend here its application. We suggest that instead 
of deforming the second term in (|1.3| ), by the same principle one may equally well deform 
the last term or possibly all terms. We shall demonstrate that the former case possesses 
some advantageous features when compared with the previously outlined deformation. 
Let us start by using the same PT-symmetric deformation principle 

u x -> -i{iu x f e G E (2.1) 

as employed in 0, albeit now for the last term. This amounts to replacing the third 
derivative as 

u xxx -» ie{iu x f~ 2 [(e - l)u 2 xx + u x u xxx ] . (2.2) 
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In this way, simply applying ( |2.2[) to flLj), we obtain a new PT-symmetric deformation of 
the KdV-equation 

u t + uu x + ze(e - l)(iu x ) £ ~ 2 u 2 xx + e(m a .) e "" 1 « a!xa! = 0. (2.3) 

At first sight the deformation ( |2.3| ) appears to be far less appealing than the deformation 
Ql.4| ). In the latter the effect of the deformation was simply that the nonlinear term of the 
KdV-equation has become somewhat more nonlinear, whereas in ( |2.3D we have replaced 
the linear term by two highly nonlinear terms. Nonetheless, as a trade off the deformation 
Q2.3] ) has some very attractive features, which are not present in (|F4|). For instance, having 
a physical application in mind we expect the deformed equation to be at least Galilean 



invariant just like its undeformed counterpart ( |l.2[) . This property is lost in (1.4), but 



instead (|2.3| ) is Galilean invariant, as it remains invariant under the transformation 

x — > x — ct, t — > t, u — > u + c, (2-4) 
where c is the velocity of the moving reference frame. Furthermore, it was difficult to 



construct conserved quantities for (1.4). Only two charges could be constructed so far and 
in addition they turned out to be complicated infinite series. We shall now demonstrate that 
this task is surprisingly simple for ([2.3|), despite its high degree of nonlinearity by relating 
it to a Hamiltonian formulation, which seems also impossible for ( |1.4| ) as it appears to be 
a non-Hamiltonian dynamical system. 

3. PT-symmetric deformations from a Hamiltonian formalism 



As we remarked, the PT-symmetry analysis, which led to (1.4), was carried out directly for 



the equation of motion. Recalling that (1.1) was obtained as a deformation of the standard 
harmonic oscillator and that this principle has been applied to various Hamiltonian systems, 
it appears highly desirable to perform deformations for the KdV system also on the level of 
a Hamiltonian. This will enable us to relate these systems to the arguments, which allow 
statements about the reality of the spectrum by utilizing "PT-symmetry as outlined in the 
introduction. In the equation of motion this property enters more indirectly and it is less 
clear which kind of conclusions can be drawn from the symmetry property. 

It is well known for a long time that the KdV-equation can be formulated as a Hamil- 



tonian system [ |14| , |i5| , 16, 17]. Thus in this spirit and in more direct analogy to the 



construction of Ql.l| ), we propose to study the new non-Hermitian Hamiltonian density 

H =u 3 - j^( iu xY +1 £ G M. (3.1) 

For £->lwe recover the standard Hamiltonian density for the KdV-equation. Clearly 7i 
in ( |3.1| ) is PT-symmetric, since it remains invariant under the transformation: t —* —t, x — ► 
—x,i — » —i and u — > u. Similarly as in the standard quantum mechanical setting, outlined 
in the introduction, PT-symmetry can be utilized to ensure the reality of the energy E, 
which follows trivially with TL{u{x)) = Tv(u{— x)) 

E= H(u{x))dx = - I H(u{-x))dx= / 7i\u(x))dx = E^ (3.2) 
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Let us now derive the corresponding equation of motion by invoking the variational 
principle for the Hamiltonian H(u) = J Tidx 



du _ d (8H{u)\ d fSj Hdx\ _ d /^°o d n dU 



\ d ( ■^—•.oo 

) = dx \^n=0^ J 



(3.3) 



dt dx \ 8u J dx \ 8u J dx \ 'n=o dx n du r . 
Evaluating (fTf) for H in yields 

u t + (-3-u 2 + eiiu^- 1 u xx ) x = 0, (3.4) 

or when not written MS M conservation law 

u t — 6uu x + ie(e — l)(iu x ) £ ~ 2 u xx + £{iu x ) e ~ l u xxx — k = 0, (3-5) 



where k is a constant. Note that ( |3.5| ) is almost ( p.3[ ), but corresponds to a deformation of 
the scaled KdV equation (|1.3j), with a = /3 = 1, 7 = —6 and k = 0, which, depending on 
the context, is also frequently used in the literature for convenience. 

3.1 Integrals of motion and conserved quantities 

Having seen how to obtain the PT-symmetrically deformed KdV equation (|2.3|), or more 
precisely its scaled version (|3.5|), from a Hamiltonian principle, we shall demonstrate next 
that is has further interesting properties, which are absent in the deformation ( |1.4| ). As 
mentioned, for (1.4) the authors of could only construct the two first conserved quantities 
in form of complicated infinite sums. Here we find instead that for ([II]) these quantities 
can be computed in a straightforward manner. Assuming to have a conserved quantity of 
the form Z^ n > = J T^ n 'dx, all we have to verify is whether its Poisson bracket with the 
Hamiltonian is vanishing, see e.g. [jjlf l. Viewing T^ n \u) and H(u) as functionals of u we 
have by definition 

Let us now employ Q3,6|) to establish that 



= J u dx, J (2) = J u 2 dx and T (3) = H(u), (3.7) 

are indeed preserved under an evolution in time. We find that these quantities are conserved 
when we invoke as standard boundary condition the non-compact or compact case for 
u,u x , . . ., that is being either vanishing at infinity or periodic in space, respectively. This 
is easily seen by computing 

— = if} = / (3n 2 - eiiu^- 1 u xx ) x dx = 0, (3.8) 

{j(2), if} = J ^ 3 - J£-(i^) £+1 - 2eu{iu x f- 1 u xx J dx = 0, (3.9) 



dt 

^ = {l^,H} = -{H,l^}=0. (3.10) 
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The last conservation law follows trivially from the anti-symmetry property of the Poisson 
brackets. We can also be more explicit and compute the corresponding flux. Constructing 
vanishing Poisson bracket amounts to seeking solutions of the conservation law 

T (n) + x (n) = Q) (3 _ n) 

with -X^ being the nth flux and being the nth conserved density. Then = 
J T^ n 'dx is a conserved charge provided the appropriate boundary conditions hold. The 
case n = 1 corresponds to the equation of motion itself as can be read off directly from 
( |3.4| ). For the case n = 2 we may re- write ( |3.9[ ) as a conservation law in the form 



{u 2 ) t + (j^- £ {iux) e+l + 1eu(iu x f- x u xx - 4u 3 ^ = 0. (3.12) 

We can also be more concrete about = Ti and compute the associated flux 

e 2 9 
* (3) = (y " z)(iu x ) 2£ - 2 u 2 xx + 3 (euu xx - 2u%) u{iu x ) £ - 1 - ie{iu x f £ - 1 u xxx - -u\ (3.13) 

thus confirming ( |3.10[ ). At this stage it is not clear whether there exist higher conserved 
quantities. However, we suspect that similarly as for most cases of the modified KdV 
equations and the generalized KdV equations only three charges exist. We recall that the 
equation u t +u p u x +u qx = is only integrable, i.e. possesses an infinite amount of conserved 
quantities, for the cases q = 3,p = 1, 2; g = l,p£N and g£N,p = l, see e.g. ||. 

3.2 Solutions of the equations of motion 

We shall now construct solutions of the equations of motion Q3.5| ). One may expect to find 
a rich variety of different types of solutions similarly as for the standard KdV equation. 
Over the years several methods have been developed to find such solutions ranging from 
minimizing the sum of the conserved charges [18|, the inverse scattering method [19]. 



Hirota's bilinearization method |2C[ |, etc. Some methods demand as a prerequisite the 
model to be integrable. As this feature is not guaranteed for the model at hand, in fact the 
conjecture is that the model is not integrable, our aim is here just to obtain a first impression 
in order to indicate that the above family of equations deserve further attention. Following 
a simple procedure which has turned out to be useful for the standard KdV equation, we 
may integrate ( |3.5[ ) directly by assuming the solution to be a steady progressing wave 

u(x,t) = w(kx — cot) = v(x — ct), (3-14) 



with c = u/k. Substituting ( p.l4| ) into the equation of motion (3J3) yields after some 
straightforward manipulations 



iTv(4n+3e+l) 

VX" J = e 



,W ,. 2(1 + 5) 



( V + 7y V +KV + K) 



i 



£ + 1 

(3.15) 



with k being an additional constant of integration and n labeling the various branches of 
the function. Separating variables then yields 

I7r( 4n + E -1) / g \ 5+1 f ([ v 

*- d = e 2 " + - (f+t) SwTW^TWn^- (3 ' 16) 
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Apart from computing the integral in fl3,16|) , the main problem is here that we need to 
solve the equation for v in order to obtain v(x — ct). This is only possible in very few 
exceptional cases, but the knowledge of the inverse function (x — ct)(v) in some domain 
will be valuable as it provides the information about the kind of general behaviour which is 
possible. For convenience we choose the dispersion relation and the constants of integration 
to be parameterized as 



c = 4k 2 (2 



m) 



4/c 4 (l — m) and k = 0. 



(3.17) 



This choice is guided by the known solutions for e = 1 and leads naturally to three quali- 
tatively different cases. 

3.3 Analogues of the cnoidal solution 

Let us first recall how to solve equation ( |3,16| ) for the case e = 1, which should result into 
an elliptic integral as we integrate the inverse of the square root of a cubic polynomial. 
With the choice of constants ( 3.17 ) we may bring ( |3.16 ) into the usual form of an elliptic 
integral 



kx — ut = ±—= 
V2 



dt 



-2k 2 



y/t 3 + 2k 2 (2 - m)t 2 + 4fe 4 (l - m)t 



4>{w) 

I 7T 



de 



(3.18) 



msm 



with (p(w) = arcsin t/(1 + w/2k 2 )/m. From (|3.18|) we deduce therefore that w{kx — ct) 
becomes the well known cnoidal solution for the KdV equation 



u(x, t) = —2k 2 dn 2 (kx — uot\m), 



(3.19) 



with dn being a Jacobian elliptic function depending on the parameter m G [0,1], see 
e.g. for notation and properties. As ( 3.19[ ) indicates for e = 1, the cases m = 0, 1 are 
special in general. For generic values of e we evaluate ( |3.16 ) with the parameterization 
d3l7D to 



x — ct 



i7T (4n + e — 1) 

e 2 (!+ £ ) 



u(l + e) 



1 



1 



1 



4fc 2 (l - m) 
l + 2e -v 



i 

l+e 



(3.20) 



,1 + E'l + e'l + e' l + e ' 2/c 2 (l-m)' 2k 2 ) ' 
Here F\ is the Appell hypergeometric function defined via a double infinite sum as 



Fi(a;(3,/3';~f;x;y) 



(«)n+m (P)n (P')m ^n^m 



m=0 n=0 



n!m!(7) 



(3.21) 



n+rn 



with (a) n := Hfc=i( a + k— 1). Since we can not solve ( 3.20j ) for v let us plot (x — ct) as a 
function of v and search for real solutions. 

We depict our findings in figure 1. For e = 1 we recognize the cnoidal solution ( 3. 1S| ) ■ 
For clarity we did not indicate the vanishing imaginary part in this case. For the other 
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values of e we find always two different types of solutions. The first resembles qualitatively 
the cnoidal solution and is either real for v G [— l/2/c 2 ,0] or v G [0, l/2k 2 ]. In figure 1 we 
present e = 3; n = 2, 4; k = l/y/2 for the former case and e = 5; n = 2, 5; k = i/v2 for the 
latter. The second type is more similar to the tan 2 solution for e = I to be discussed in 
the next section. These solution are real either for v G (— oo,0] or v G [0, oo). In figure 1 
the former case is illustrated by e = 3; n = 2, 4; k = i/y/2 and the latter bye = 5; n = 2, 5; 
k = l/y/2. 




1 1 1 1 1 1 1 1 1 1 
-~ / 


6= l;n= 1,2 

c = 3;n = 2,4 

e = 5; n = 2,5 

e = 3;n = 2,4 
s = 5; n = 2,5 


s . ~~~ 

i . i . i i . i . i i 



1.00 -0.75 -0.50 -0.25 



0.00 
V 



0.25 0.50 0.75 1.00 



Figure 1: (x-ct) as a function of v for e = 1, 3, 5 for some particular branches. 
3.4 Analogues of the tan 2 solution 

Next we consider the limit m — > 0. Keeping the parameterization (|3.17|) we make the 
further convenient choice k = dbl/v2, similarly as in the previous section, which amounts 
now to the boundary condition k = 1. Then the Appell hyper geometric function F\ reduces 
to the Gauss hypergeometric function 2-Fi defined as 



iF\ (a; P; 7; x) := ^ — — — x 



n=0 



n!(7) n 



Fi(a;/?/2,/?/2; 7 ;z;x) 



Using furthermore the identity 

2F1 (a; 2/3; 2a + /3; s) = ax~ a B x (a, 1 - 2/3) , 
where (a, /?) is the incomplete beta function 

B z (a,P)= f t a ~\l-tf-\ 
Jo 

we obtain the simpler expression 

1 

e e — 1 



(3.22) 



(3.23) 



(3.24) 



?7r(4n + £ — 1) 

x -4t = e 2 ( 1 + £ ) 



e + 1 



1 



e+l'e + 1 



(3.25) 



For e = 1 ( 3.25 ) reduces further to x — 4t = \/2 arctan(±- v /u), which may be solved 
for v, such that we obtain u(x,t) = tan 2 [(x — At)/y/2\ as a solution for the standard KdV 
equation. For generic values of e we depict (jjT25p for various values of the parameters 
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in figure 2. For e = 1 we perceive the real solution x — 4t = y/2 arctan(±y / w) in panel 
(a). A qualitatively similar type of solution is obtained for instance for some branches 
for e = 5 as is seen also in panel (a). Panel (b) confirms that for v > this solution 
is real. (The solid line is on top of the dashed line) Interesting qualitatively different 
types of solutions are obtained for instance for some branches for e = 3,11. We observe 
from panel (c) that these solutions are very reminiscent of the one soliton solution, to be 
discussed in the next section, albeit with the fundamental difference that they are not 
vanishing asymptotically for large (x — ct). This is seen simply by using the property 
B\ (a,f3) = T(a)T(P)/T(a + (3) of the incomplete beta function. For v-tlwe obtain in 

. . „r(4n + e-l) / \ 1/0 + 1) 1 „ , 

( PT25D the definite values e 2(i+ e) ^ r ( ^ ) r ( |^ ) /T ( ^ ) . This limit is finite 

for the parameter range except for e = 1, when lim x ^,oT(x) — > oo. In this case we obtain 
a purely complex one soliton solution as can also be seen clearly in panel (b). Having 
Galilean invariance for our equations, we may also move this function as v — ► v + 1, such 
that the tails are located at v = rather than v = — 1, which is a more familiar setting. 

3 
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Figure 2: (x-ct) as a function of v for e — 1, 3, 5,11 for some particular branches and m=0. 



3.5 Analogues of the one soliton solution 

Next we take the limit m — * 1 corresponding to the special case k = 0, which implements 
vanishing boundary conditions. Indeed, adopting the parameterization ( 0.17 ) the limit 
m — ► 1 in fl3.19|) for e = 1 yields the asymptotically vanishing single soliton solution 
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u(x,t) = — 2A; 2 sech 2 (Ax — ut). Taking this limit in ( 3.20| ) for generic values of e gives 



i7r(4n + £ — 1) 

x-ct = e 2 ( 1+£ > 



e + 1 



(2k z )~ B_ 



\e + l e + 1 



(3.26) 



We depict this function for various values of the parameters in figure 3. The famous 
one soliton solution is clearly visible for e = 1. In the other cases we obtain again two 
qualitatively different types of real solutions. One type being real in the finite ranges 
v G [-1/2A; 2 ,0] and v G [0, 1/2A 2 ] exemplified by e = 5; n = 1,4; A; = l/y/2 and e = 3; 
n = 2,4; A; = i/V2, respectively. The other type is real in the ranges for v G (— oo,0] 
and v G [0, oo), which we illustrated in figure 3 by e = 3; n = 1, 3; A; = i/\/2 and e = 5; 
n = 2,5; A; = l/y/2, respectively. 
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Figure 3: (x-ct) as a function of v for e = 1, 5 and e = 3, 11 with fe = ±l/y/2 and fe = ±i/V%, 
respectively, for some particular branches with m=l. 



4. Conclusions 

Alternatively to J/]], we proposed a new "PT-symmetric complex deformed version of the 
KdV equation. The suggested deformation allows for a simple non-Hermitian Hamilto- 
nian formulation involving a Hamiltonian density very reminiscent to the prototype VT- 



symmetrically complex deformed quantum mechanical system (hi). The model ( |3.1[) is 



Galilean invariant and three charges, related to the conservation of mass, momentum and 
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energy, together with their conservation laws, were constructed. We demonstrated that 
there exist steady progressing wave solutions for these models and identified analogues to 
the cnoidal and tan 2 solution. However, we did not find asymptotically vanishing analogues 
to the one soliton solution. 

Clearly there are many important questions left to be answered. It would be interesting 
to establish that there exist three and only three charges for the proposed deformation. 
Besides solving the equations more explicitly it will be natural to seek for solutions on 
some rays in the complex plane. It will be straightforward to extend these considerations 
to the modified KdV and the generalized KdV equations. We shall leave these issues for 



future investigations [22 |. 
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